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We derive a formula, useful for first-principles calculations, 
which relates the free energy of an oxide/metal interface to 
the free energies of surfaces and the work of separation of 
the interface. We distinguish the latter mechanical quantity 
from the thermodynamic work of adhesion, and we describe 
explicitly how both may be calculated. Our formulae for in- 
terfacial and surface energies are cast in terms of quantities 
which can be calculated or looked up in tables, and include as 
additional parameters the ambient temperature and partial 
pressure of oxygen Po 2 - From total energy calculations for 
the Nb(lll)/a-Al 2 3 (0001) interface, free Nb and A1 2 3 
surfaces, we obtain firstly numerical estimates of the works 
of separation, which are independent of Po 2 ■ We then obtain 
surface energies, interfacial energies and the equilibrium work 
of adhesion as a function of Po 2 ■ 

PACS numbers: 68.35. G, 71.15.N, 73.20 



I. INTRODUCTION 

Oxide-metal interfaces continue to be studied inten- 
sively because of the many ways in which they are of com- 
mercial and scientific importance. Applications range 
from the nanoscale in microelectronics packaging to the 
macroscale engineering of thermal barrier coatings or the 
formation of protective scales. The science of these in- 
terfaces has been addressed in volumes of conference pa- 
pers and reviews There are also reviews in the 
literature |5||| which specifically address the theoretical 
questions about the nature of the bonding at these in- 
terfaces, such as: what determines the site-preference of 
metal atoms on the oxide surface; whether the bonding 
can be thought of as predominantly covalent or metallic 
and how to quantify these concepts; whether a simple 
classical image model can be used to interpret the bond- 
ing; what is the strength of adhesion of metal to oxide. 
The basis for answering these questions is to have reliable 
calculations of the electronic structure and total energy 
of particular surfaces and interfaces. Such calculations 
came of age over the past ten years or so with the use 
of first-principles methods. These mainly apply density 
functional theory (DFT) and the local density approx- 
imation (LDA) which are the basis of the calcu- 
lations we shall report here. Hartree-Fock calculations 
are also feasible and have been applied effectively to the 
Ag/MgO interface M, although they tend to be more ex- 



pensive than DFT for larger systems. Since the reviews 
cited, there have been numerous applications of DFT to 
study bonding in the initial stages of deposition of metal 
on oxide, with cluster or multilayer geometries, notably 
on MgO [jl0|-^l| , but to a lesser extent on more complex 
oxides such as Ti0 2 p|, MgAl 2 4 pi and a-Al 2 3 



The bonding of Nb to a— AI2O3 has long been a sub- 
ject for experimental work, because besides its relevance 
to electronic components it offers practical advantages for 
sample preparation: the two materials bond strongl y |2q | 
(anomalously strongly according to a recent study ||29| ), 
do not react chemically and have similar coefficients of 
thermal expansion. In the orientation Nb(lll)/Al203 
(0001) there is a lattice mismatch of < 2%, allowing the 
preparation of a nearly coherent interface (using molec- 
ular beam epitaxy), the atomic structure of which has 
been studied by high resolution transmission electron mi- 
croscopy (HRTEM) and analysed in detail fc]-§4|. This 
interface was the subject of first-principles calculations 
which used periodic boundary conditions, making the 
reasonable assumption that the effect of misfit disloca- 
tion can be neglected |p5|-p7|] . Our recent work [^6| anal- 
ysed the nature of the bonding in detail by calculating 
Mulliken populations and bond orders, concluding that 
the bonding across the interface is strongly ionic. The 
work of separation W sep of the interface was calculated, 
and found to be very high: of order 10 Jm -2 when nio- 
bium was bonded to the oxygen terminated AI2O3 sur- 
face. Lower energy pathways for the cleavage of this in- 
terface would be within the Nb metal or the oxide itself. 
Two other interfaces were studied corresponding to the 
two other possible terminations of bulk AI2O3 (0001); 
namely the stoichiometric, aluminium termination (one 
layer of aluminium) and the aluminium-rich termination 
(two layers of aluminium). HRTEM could not distinguish 
between the stoichiometric termination and the oxygen 
termination, however evidence from electron energy loss 
spectroscopy |35|| favoured the oxygen termination. 

We point out here that as far as we know the question 
of which termination is more stable has not yet been ad- 
dressed in all the theoretical work which has been pub- 
lished so far on any oxide-metal interfaces. The struc- 
tural predictions have been confined to the question of 
the relative displacement of the crystals, parallel and per- 
pendicular to the interface, and the local relaxations of 
atoms at the interface, as well as the energy needed to 
separate the crystals W sep . This has been done for in- 
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terfaces with different terminations or local stoichiome- 
try; all calculations were carried out with atoms at rest 
(T = OK) and minima in the total energy were located 
as a function of atomic positions. However, the question 
as to whether the oxygen terminated or the aluminium 
terminated interface is more stable was not discussed. It 
is well known that this question can only be answered 
with respect to the chemical potentials of the species in 
the environment with which the interfaces are in equi- 
librium, which is normally characterised by temperature 
and partial pressure of oxygen |3(|, and the difficulty of 
relating the quantities accessible to a first principles cal- 
culation to these parameters may have been a reason for 
leaving this question to one side. 

The main purpose of our present paper is to show how 
in fact we are already able to make predictions of the sta- 
bility of different interfaces when they differ not only in 
structure but also in composition. With certain simplify- 
ing assumptions we show how this can now be done with 
little more effort than the calculations which need to be 
done to calculate the work of separation, and we present 
first results for the Nb/A^Oa interface. The ingredients 
of the theory are the works of adhesion and surface ener- 
gies. For these we draw upon the results reported briefly 
in p6| , supplemented by some further calculations to dis- 
cuss the case of oxygen on the Nb surface. The basic 
theory is outlined in Section II. We derive the equations 
for a general A rn O n oxide in contact with a metal B\ 
it would be a short step to generalise them still further 
to an interface between arbitrary compounds. Essentially 
the same thermodynamics was applied by Wang et al pTj] 
in calculations of the surface energy of oxides with dif- 
ferent terminations, over a range of chemical potentials 
of oxygen; our theory makes the further connection to 
the temperature and in particular the pressure of oxy- 
gen, which are the parameters directly under the control 
of the experimentalist. A detailed study of the AI2O3 
surface will be reported elsewhere |3£| . 

The plan of the rest of the paper is as follows. Sections 
III-V cover aspects of our total energy calculations which 
were not dealt with, or dealt with only briefly, in our Let- 
ter ||q] . In Section III our method of total energy calcula- 
tion is summarised. In Section IV we describe the atomic 
relaxations parallel to the interfaces which are generally 
not commented upon. Although this structural aspect 
is not central to the thrust of our paper, it turned out 
that lateral relaxations also have a strong part to play 
in determining the interplanar relaxations and energies 
reported previously, and we therefore describe them for 
completeness. In Section V we describe and comment on 
the results for the work of separation on different planes 
and with different terminations of the interface. Our cal- 
culated interfacial free energies are presented in Section 
VI and we conclude in Section VII. 



II. PRINCIPLES OF CALCULATING 
INTERFACIAL AND SURFACE ENERGIES 

Let us consider the interface between metal B and an 
oxide of metal A in equilibrium at temperature and pres- 
sure (T,P). The stoichiometric composition of the A 
oxide is A m O n . We obtain the definitions of interfacial 
quantities by referring to the contents of a periodically re- 
peated supercell of area S parallel to the interfaces which 
it may contain. All extensive thermodynamic quantitites 
in the following will refer to the contents of such a super- 
cell. The interfacial energy per unit area, counting the 
two interfaces within each supercell, is given by ]39| ] : 

lint = (G mt (T, P) - N At x A (T, P) - N ot io(Ti P) 

-N Bf i B (T,P))/2S, (1) 

where Gi„t is the Gibbs energy of the contents of a su- 
percell containing two interfaces, fj, A , Mb and \xo are the 
chemical potentials of the three components, and N A , 
Nb and No are the numbers of atoms of the three com- 
ponents within the supercell. The denominator 25* oc- 
curs because there are two interfaces in the supercell, as 
required by periodic boundary conditions. Chemical po- 
tentials here are per atom rather than per mole, which 
would be the usual convention for macroscopic thermo- 
dynamics. Special cases of Eqn. @) are when either the 
metal B or the oxide is absent from the supercell, in 
which cases we recover expressions for the surface ener- 
gies of the oxide 7,40 and the metal jb respectively: 

Iao = (G S ao(T,P) - N Af i A (T,P) - N no(T,P))/2S, 

(2) 

IB = (G SB (T,P) - NbVb(T,P))/2S. (3) 

The quantities Gs A o and Gsb are the Gibbs energies 
of slabs of oxide and of metal, with free surfaces separated 
in their respective supercells by an adequate thickness of 
vacuum. We have assumed in Eqn.(||) that the metal 
surface is clean; this will suffice for a calculation of the 
energy of the interface discussed below. However, we can 
easily consider for example adding a monolayer of oxy- 
gen to the metal surface in the calculation of Gsb(T, P). 
Since there is no separate oxide phase, the number of oxy- 
gen atoms in the system, No, now resides on the metal 
surface. The contribution —NoHo(T, P) must be sub- 
tracted in the calculation of the corresponding surface 
energy just as in the calculation of the interface energy. 
We shall in fact make this calculation of an 'oxidised' Nb 
surface in the course of obtaining the work of separation 
of an interface by a pathway which leaves oxygen on the 
exposed metal surface. 

The motivation for calculating "fi nt is as follows. An in- 
terface between two crystals requires five parameters for 
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its macroscopic specification, for example three to specify 
the relative crystallographic orientation of the materials 
and two more to specify the orientation of the interface. 
We note in passing that a free surface in contact with 
vapour or liquid only requires two parameters to specify 
its crystallographic orientation. There is always a large 
set of hypothetical interfaces which have the same five 
macroscopic parameters but which differ in their atomic 
structure and local stoichiometry. The member of this 
set which minimises ji n t for given chemical potentials 
is the equilibrium interface. So provided we know the 
chemical potentials of the components, we could in prin- 
ciple predict the atomistic structure of the equilibrium 
interface, including its local stoichiometry, by evaluating 
7i„i for each member of the set. In practice, of course, 
we can only calculate jmt for a small subset of the entire 
set, and rely on our experience and intuition, together 
with experimental information, to ensure that we have 
not omitted an important structure. Prior to the present 
work we and others have calculated total energies for a 
number of structures at OK, in which the atomic positions 
are relaxed by energy minimisation. The equations to be 
derived below show how to go the two important steps 
further, namely to correct the OK information to finite 
temperature and to take account of local stoichiometry. 

The relationship between local stoichiometry and ji nt 
is well known in thermodynamics as the Gibbs adsorp- 
tion equation, in which local stoichiometry is measured 
in terms of excesses Ti of one or more components la- 
belled i. Our final version of Eqn. (|l|) will be in terms of 
the excess of oxygen at the interface with respect to the 
metal A, per unit surface area, which is defined as: 

ft 

To = (No - ~N A )/2S. (4) 
m 

This choice of component i is arbitrary; we could 
equally well work in terms of Tai, because they are re- 
lated through 

mTo + nT A i = 0. (5) 

We note three further points in connection with ex- 
cesses. Firstly, each further component in the system 
would introduce another excess, each excess being re- 
ferred to the same designated component. Secondly, a 
stoichiometric interface is by definition one for which all 
the excesses vanish. Finally, since one of the phases is the 
pure metal B, there can be no excess of the metal B. In 
particular, if B were in fact also Al, thereby reducing the 
number of components to two, the interface could not be 
described as having an excess of O or Al. For more dis- 
cussion of the thermodynamics of excess quantities the 
reader is referred to |39] . 

A difficulty up to now has been to calculate the chem- 
ical potentials involved in these equations and, more 
specifically, to relate them to given experimental con- 
ditions. In the following we show how Eqn. (Ill) can be 



reformulated to relate ji n t to the partial pressure of oxy- 
gen P 02 . 

First we define the work of separation W sep of the in- 
terface. It does not refer to chemical equilibrium states 
and therefore does not involve chemical potentials of the 
separate components: 

W sep = (Gsao + Gsb — G int )/2S 

= 1AO + IB - lint ■ (6) 

For brevity we do not explicitly indicate the temper- 
ature and pressure dependence of all the terms unless 
it needs to be emphasised. An important point to note 
about this quantity, which makes it relatively straightfor- 
ward to calculate, is that the separate slabs of metal and 
oxide have exactly the same composition as the two slabs 
which are joined to form an interface. This would not in 
general be the case if these slabs and the interface were 
in equilibrium with a given environment (constant /Ltj), 
because one would expect for example some loss or gain 
of surface oxygen or metal from the oxide to the vapour 
phase when the surfaces are created. If the interface as 
well as the exposed surfaces are the ones which are in 
chemical equilibrium (which we denote by superfix eq), 
an equation similar to (|^) defines the work of adhesion: 

W ad ^i e A q +lB q -7Z- (7) 

This is the quantity of relevance to contact angles and 
wetting for example, and unlike W sep it is not obtainable 
by a simple comparison of three total energies. 

Calculations of W sep for Nb/A^Oa were reported in 
our Letter |2{| and these have been extended here, as 
described in the following sections. W sep is probably 
more relevant than W a d in formulating a fracture cri- 
terion, when internal surfaces are formed which are not 
in equilibrium, but in order to predict the equilibrium 
structure of interfaces we also need to be able to evalu- 
ate Eqns.(|)-@. 

We now introduce the quantity gAO-, the Gibbs energy 
per formula unit of bulk A m O n in equilibrium with metal 
A and oxygen in vapour form: 

g AO [T, P) = TnfiA (T, P) + nfi (T, P) , (8) 

so that 

Gao = {N A /m)g A o (9) 

is the Gibbs energy of a stoichiometric cell containing 
Na atoms of A. Inserting (@) and (0) into (||) gives the 
surface energy of the oxide in a form which makes the 
effect of the excess oxygen explicit: 

Iao = (G S ao~Gao)/2S-To^o- (10) 

Consider now how to go about calculating the two sur- 
face energies from (§) and @, which we will eventually 
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combine with W sep in (||) to give us ji nt . The Gibbs en- 
ergy of all slabs can be calculated at T = OK and P = 
from first principles, it is just the total energy. If the 
slabs are bulk pure material, their Gibbs energy can be 
corrected to temperature T by using experimental spe- 
cific heat data. On the other hand when the slabs are 
separated in the supercell by a layer of vacuum to repre- 
sent free surfaces, there is no such experimental data and 
the correction to finite T could be done by calculating the 
phonon spectrum and using the quasiharmonic approxi- 
mation for the free energy. This has been done previously 
for classical ionic models by Taylor and coworkers p0| , 
in order to obtain the temperature dependence of their 
surface energy, but we have not yet made the equiva- 
lent calculation with our ab initio code. For a metal slab 
(Ag), the quasiharmonic free energy based on ab initio 
phonon frequencies was recently calculated by Xie and 
coworkers EJ. In the case of the pure metal slab, the 
chemical potential /is is the Gibbs energy per atom of a 
bulk slab. The surface energy of B is therefore obtained 
from the results of two supercell total energy calculations 
in the standard way. The main present issue, which is 
less familiar in the context of total energy calculations, 
is how to calculate the significant term due to the chem- 
ical potential of oxygen, which must be included when 
the surface of the oxide is non-stoichiometric (To 7^ 0) . 

The chemical potential of oxygen is well described in 
terms of its partial pressure Po 2 by the standard ideal 
gas expression 



& + \kt log(PojP°) 



(11) 



In Eqn. ([ll]), /1q is the oxygen chemical potential in its 
standard state (STP) at T°=298.15K, P°=lat. Chemists 
would set /iQ to zero by definition, but we cannot do that 
since our zero of energy is already defined as the energy 
of separated ions and electrons at T = OK. On the other 
hand the energy of oxygen molecules is not something 
we want to calculate, since there are well known prob- 
lems in using density functional theory for this system. 
Fortunately, we can circumvent the problem by using a 
thermodynamic cycle. From the defining equation for the 
standard Gibbs energy of formation AG^ : 



9°AO=mfJp A + nfi o + AG AO , 



(12) 



we obtain the troublesome oxygen chemical potential at 
STP in terms of g AO and AG ao . The quantities g AO 
and /j, a are things we can calculate accurately, and we 
can look up AG AO in tables of thermodynamic data. 

Inserting fi% from © into (0) and @ into (0) 
gives us our final expression for the surface energy of the 
oxide: 

Iao = (G S ao(T,P) - ^g A o(T,P))/2S 
m 

-To(g AO ~^A~^G AO )/n 



-To-kT log(Po 2 /P°). 



(13) 



from which we obtain the final expression for the inter- 
facial energy by subsituting (|l^) into (||): 



lint =J B {T,P)-W sep {T,P) 



<G SAO (T,P) 



N A 



9ao(T,P))/2S 



-T {9 A0 -mii A -AG A0 )/n 
-Vo\kT \og{P jP ). 



(14) 



The quantities g AO and /i A entering the third line of (|TJ) 
are well described by T = OK quantities which we calcu- 
late. It can be verified that correcting them to standard 
state has a negligible effect on the surface energy. 

The minimum physically meaningful value of Po 2 , 
which we denote Pq™, is set by the condition that if 
Po 2 < Po" 1 t ne oxide would spontaneously decompose 
into metal and oxygen. Neglecting the small variation 
in solid energies with temperature by comparison with 
AG AO this condition is: 



log(P™"/P°) = 



nkT 



-AG 



AO- 



(15) 



Similarly, the maximum physically meaningful value of 
Po 2 is defined by the lowest standard Gibbs energy of 



formation of a metal B oxide AGg . 



iog(P«r7 p0 ) 



'kT 



AG" 



BO 



(16) 



where the first oxide to form would have the stoichiom- 
etry B m 'O n >. The thermodynamic data used here are 
summarised in Table || 



III. METHOD OF TOTAL ENERGY 
CALCULATION 

For the interface calculations we use the total energy 
plane wave pseudopotential method based on Lanczos di- 
agonalization of the Kohn-Sham density matrix pffij . The 
supercell has the form of a rhombohedral prism and in 
the stoichiometric slab it contains 45 atoms: 14 Al, 21 O 
and 10 Nb atoms (see Fig.[l]). By stripping off the outer 
plane of Al from each interface we obtain an interface 
which is O-terminated with an O excess To • S = +1.5 
atoms per surface unit cell. By adding the surface plane 
of Al atoms to the neutral interface we make an oxygen 
poor interface, with the negative O excess To ■ S = —1.5. 
The total energy of the contents of a supercell is min- 
imized with respect to the atomic coordinates by the 
quasi-Newton method with Hessian updated using the 
Broyden-Fletcher-Goldfarb-Shano (BFGS) method. The 
pseudopotential for Nb was of Troullier-Martins form 
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p4j , with s and d non- locality. The pseudopotential for 
O was of optimised form |0, with s non- locality. The 
pseudopotential for Al was of Gonze type with s non- 
locality. 

All calculations were made with two ^-points in the ir- 
reducible wedge of the Brillouin zone, and with a plane- 
wave cut-off of 40 Ry. The effect of increasing the plane 
wave cutoff from 40 to 60Ry was to reduce W sep by 3.3% 
for the Nb/Al interface, which we take as a satisfactory 
indication of the basis set convergence. For the neutral 
45 atom interface we have made test calculations with 
six fc-points which results in a decrease of total energy 
by about 4mRy and very small (< 10 _3 nm ) changes 
of relaxed positions of atoms compared with two /c-point 
calculations. The effect of increasing the /c-point sam- 
pling from 2 points to 9 is to change W sep by less than 
1%. 

By doubling the original unit cell in the x-y plane we 
obtained a 180 atom cell, with which we recalculated the 
wavefunctions at the gamma point with the previously 
relaxed atomic coordinates. The gamma point wavefunc- 
tions in this cell were used for Mulliken population anal- 
ysis which was made by projecting the optimized wave 
functions onto the pseudoatomic orbitals \(f>i a > (i labels 
site, a - orbitals) according to the procedure suggested in 
p7| . The "spillage" of each occupied orbital ip was less 
than 1.5%. 



IV. RELAXATION OF THE INTERFACE 

The slab with which the Nb(lll)/Al 2 O 3 (0001) inter- 
face was modelled is shown in Figjl]. The interlayer relax- 
ation of the interface has been reported previously [^6| , 
and we refer to that paper for results. Here we mention a 
feature which has not previously been discussed, namely 
the relaxations parallel to the interface, which we refer to 
as in-plane relaxations. It has been found that to make 
a calculation of the interlayer relaxation of the alumina 
surface one needs to take into account the in-plane re- 
laxations of the oxygen atoms |^(| , which were neglected 
in some earlier work [^5| . The present results show that 
in-plane relaxation of the oxygen ions is a general fea- 
ture of the structure near the interface. The geometry 
of these relaxations is described by two parameters; the 
rotation a and the bond length extension 8r of the equi- 
lateral triangle of oxygen atoms in a plane, illustrated 
in Fig.||. The rotation and dilation of this equilateral 
triangle does not break any symmetry, preserving for ex- 
ample the three fold axis about the centre of the triangle 
concerned. 

The calculated values of a and 8r are shown in Table || 
for the terminating and second layer oxygen planes (de- 
note by subscripts 1 and 2) in five cases. The first two 
cases, labelled N(b)/A(0) and N(b)/A(A1) are the O and 
Al terminated bulk Nb/Al 2 03 interfaces. Cases three 



and four are pure alumina surfaces, labelled A(O) and 
A(A1) to indicate that they are oxygen and aluminium 
terminated respectively. Case five, labelled N(m)/A(0) 
is a Nb monolayer on an oxygen terminated alumina sur- 
face. One can see that the in-plane relaxation is a feature 
of all the systems studied. From the evidence of the first 
two layers, the rotation of O-triangles and the increase 
of the 0-0 bond lengths appears to be localised near 
the surface of the alumina. The Nb monolayer on the 
O-terminated surface of alumina shows this effect most 
strongly, which is quite surprising, since the interplanar 
relaxation in this case is much less than that of the Al 
layer for which the Nb substitutes. There is experimen- 
tal confirmation of the effect, obtained by small angle 
X-ray diffraction flgf] , in the case of the A(A1) surface. 
In this case the experimental results are 8r\ — 4.5% 
and ai—6.7°, compared with our calculated results of 
Sri = 3.2% and ai=3.1°. The agreement is only quali- 
tative. 



V. WORK OF SEPARATION 



All our results for the calculated work of separation 
W sep of different interfaces and cleavage planes are shown 



in Table III. The column of 'unrelaxed' results refers 



to values obtained by assuming bulk unrelaxed atomic 
positions both at the interface and for the free surfaces. 
The interplanar spacing between O and Nb across the 
interface in this case was simply taken as the bulk spacing 
between O and Al planes. The results in the 'relaxed' 
column are calculated with atomic positions relaxed both 
before and after cleavage. 

The effect of relaxations on W sep naturally depends 
on the interface. It is most pronounced when an Al ter- 
minated AI2O3 surface is exposed, because of the large 
relaxation of this surface, which lowers its surface energy 
by about 1.5 Jm -2 . On the other hand, if the relaxation 
of the interface dominates the energy balance in Eqn.(|^) 
then the relaxed value of W sep c&n even be larger than 
the unrelaxed value, as in the case of cleavage between 
Nb and O at the N(b) /A(O) interface. 

The lowest value of W sep , 2.7 Jm~ 2 , is found for 
the cleavage of bulk Nb from the stoichiometric Al- 
terminated alumina. The highest values are found for 
the cleavage of bulk Nb from the O-terminated alumina 
surface. Indeed we can deduce from Table [II that this 
interface would be unlikely to separate between Nb and O 
planes, but would prefer to separate inside the Nb, leav- 
ing a monolayer of Nb on the surface, or even between O 
and Al, leaving a monolayer of O on the Nb surface. 

The highest value (relaxed) of all in Table I is for the 
cleavage of pure AI2O3 between O and Al planes. Exper- 
imentally, a— AI2O3 does not cleave on the basal plane 
at all, but its lowest energy cleavage on this plane would 
clearly be between Al planes. This is what one expects 
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on the basis of charge neutrality arguments, because by 
cleaving between Al planes two identical, neutral surfaces 
are created. On the other hand by cleaving between O 
and Al, different surfaces are created which, in order to 
be neutral, require the oxygen or aluminium at the sur- 
face to be in an unfavourable valence state, hence this is 
a final state of especially high energy. 

The above interpretations are supported by the Mul- 



liken populations shown in Table IV for the three inter- 
faces between bulk Nb and AI2O3. We make the usual 
caveat here that Mulliken charges do not have absolute 
significance, since they depend on the choice of basis set, 
but they are nevertheless a useful indicator of trends in 
ionicity or covalency. The interfaces in Table III are 
labelled by their oxygen excesses, to highlight certain 
trends with the stoichiometry of the interface. Bulk oxy- 
gen carries a Mulliken charge of -1, and for the oxygen 
plane nearest the interface this value is reduced to -0.99, 
-0.93 and -0.86 in turn as the excess of oxygen at the in- 
terface is increased from negative to positive. The change 
is rather modest, indicating that oxygen does not read- 
ily alter its valence state. The charge on the interfacial 
oxygen is provided by the terminating layer of Al in the 
case of Al termination, or in the case of the oxygen ter- 
minated interface the electrons are provided mainly by 
the first two layers of Nb. 



VI. INTERFACIAL FREE ENERGY AND 
OXYGEN PRESSURE 

Five surface energies are shown in Fig.[|a as a function 
of Po 2 ■ The x-axis is appropriate to a temperature of 
1500K; to obtain the results at temperature T the num- 
bers on the x-axis should be multiplied by 1500/T. The 
O-rich (T ■ S = 1.5) and O-poor (T ■ S = -1.5) alu- 
mina surfaces have negative and positive slopes respec- 
tively, while the stoichiometric AI2O3 and pure Nb(lll) 
surface energies are constant, and by chance nearly equal. 
The most negatively sloping surface energy we have plot- 
ted here refers to the Nb(lll) surface with an attached 
monolayer of oxygen. It becomes negative at an oxygen 
pressure inside the regime of stability of NbO. 

The interfacial free energy from Eqn.(|l]) is shown as 
a function of Po 2 in Fig]| c for three interfaces, O-rich, 
stoichiometric and O-poor (To ■ S = 1.5,0,-1.5). To 
generate the interfacial free energies one has to subtract 
the work of separation shown for convenience in Fig.||b, 
from the sum of the equilibrium surface energies of the 
two corresponding free surfaces. With increasing Po 2 the 
O-rich interface becomes increasingly stable, the Al-rich 
interface less stable and the free energy of the stoichio- 
metric interface remains constant, exactly parallel to the 
behaviour of the free surfaces. The interfaces can only 
be in thermodynamic equilibrium in the range of oxygen 
pressure which is indicated on the figure. At 1500K this 



range is as given in Table || at values of Po 2 above this 
range, the Nb would oxidise to NbO, and below it the 
alumina would decompose. 

The work of adhesion at a given Po 2 can also be esti- 
mated from the results on this graph using Eqn.(Q). The 
free surface and interface energies should be those with 
lowest free energy at the given oxygen pressure, and these 
can be read off from Figs.||a and||c. The result is plotted 
in Fig.gd. 



VII. DISCUSSION AND CONCLUSIONS 

We have made a careful distinction between work of 
separation, a mechanically defined quantity, and work of 
adhesion, a thermodynamic quantity, focusing on how to 
go about calculating these quantities within an atomistic 
model. We particularly consider the interface between a 
metal and an oxide, since it is of practical importance 
and since oxygen is a troublesome component for which 
to calculate the chemical potential, a key quantity in in- 
terfacial energies. A useful practical equation for the 
free energy of an interface involving oxygen has been de- 
rived, namely Eqn.(^), which gets around the previous 
difficulty by using a thermodynamic cycle to express the 
result in terms of quantities which can be readily calcu- 
lated, namely the total energies of slabs, and quantities 
which can be obtained from tables, namely the standard 
Gibbs energy of formation of the oxide and the Gibbs 
energies of the bulk materials relative to their T = OK 
values. 

To illustrate and apply the method we have made a 
number of first principles calculations for Nb(lll)/Al203 
(0001) interfaces, oxygen rich, oxygen poor and stoichio- 
metric, and for several surfaces. We fully relax the atomic 
positions in supercells using a plane wave, pseudopoten- 
tial methodology. The relaxations are significant, and in 
all cases they involve in-plane as well as interlayer re- 
laxations of the oxygen ions. Results on the work of 
separation of these interfaces were given in a Letter ptf l 
recently, and we have extended them to include the pos- 
sibility of a cleavage of the O-terminated interface which 
leaves the Nb coated with oxygen. This turns out indeed 
to be a lower energy mode of separation (4.9 Jm~ 2 ) than 
the alternative which leaves an oxygen rich AI2O3 sur- 
face behind (9.8 Jm~ 2 ), because the favourable degree of 
ionicity of oxygen is thereby preserved as it is in both 
bulk alumina and its stoichiometric surface. Consider- 
ing further the strongly bound O-terminated Nb/Al203 
interface, it turns out that the hypothetical processes of 
(i) cleavage within bulk Nb (4.2 Jm~ 2 ), or (ii) leaving 
a monolayer of Nb on the oxide surface (3.8 Jm~ 2 ), or 
even (iii) cleavage within bulk AI2O3 (3.9 Jm~ 2 ) are all 
marginally of lower energy than the cleavage which takes 
oxygen with the niobium. 
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By combining the results of our calculations with ther- 
modynamic data we obtain surface energies and inter- 
facial energies as a function of oxygen partial pressure 
and temperature. An approximation we make here is in 
omitting the temperature dependence of the solid state 
free energy, but we include the kT \og(Po 2 / P°) term 
which describes the temperature dependence of the oxy- 
gen chemical potential; this is also the term which de- 
scribes the dependence of all the interfacial and surface 
free energies on oxygen pressure. It is clear how a more 
accurate calculation could be made by implementing the 
quasiharmonic approximation to correct solid surface free 
energies, and it will probably become a routine mat- 
ter to include such a correction in future work. An- 
other approximation is made by considering only a small 
set of possible interface and surface compositions which 
we think are representative. Nevertheless, despite the 
present simplifications, some clear results have emerged. 

Of the free surfaces of AI2O3 , the stoichiometric one, 
terminated by a single layer of Al, is the most stable 
over the whole range of oxygen partial pressure up to 
over one atmosphere. It may be that a treatment of the 
temperature dependence of the energy of the slabs could 
modify the upper and lower bounds on pressure some- 
what. Correction of the LDA error is also likely to lower 
surface energies by 10-20% (I. G. Batyrev, unpublished). 
For example, work of C. E. Sims et al (49| with classi- 
cal potentials indicates that the surface energies of Al- 
terminated AI2O3 can be reduced by up to 0.2-0.3 Jm -2 
at 1500-2000K. We expect an oxygen terminated surface 
to be stable at a pressure not too far from atmospheric, 
but we cannot unfortunately be more quantitative in the 
prediction at this stage. At very low oxygen pressures it is 
also reasonable that the experimentally observed Al-rich 
(\/31 x vol) structure is stable; we cannot model a super- 
cell of the size needed to calculate this. Instead we mod- 
elled a much simpler Al-rich interface, which is predicted 
to become the most stable one just above Pq™ where 
AI2O3 decomposes. Since the experimental x V31 is 
a very Al-rich surface (To • S = —7.5 in the present nota- 
tion), the slope of its surface energy versus log(Po 2 /-P ) 
is correspondingly very steep and positive, and it must 
intersect all the other surface energies just above Pq 171 . 

The Nb free surface should obviously become unstable 
with respect to some adsorption of oxygen when P<p 2 > 
pmax^ pressure at which NbO begins to form. The 
particular configuration and concentration of an oxygen 
monolayer which we have calculated is not likely to be 
the optimum configuration of the first oxygen covered 
Nb(lll) surface, but it does become more stable than the 
free surface at pressures somewhat above Pq" x (Fig]3|a) . 

A significant new result is the theoretical analysis of 
the thermodynamic stability of the O-terminated inter- 
face, the strong bonding of which we discussed above. 
No interface is thermodynamically stable above the (very 



low) oxygen pressure at which NbO forms, but over most 
of the range below this the O-terminated interface is less 
stable than the Al-terminated one (Fig.||c), despite its 
strong bonding. In fact at the very lowest pressure of oxy- 
gen, as would pertain in the presence of pure aluminium, 
our prediction is of an Al-enriched interface. The exper- 
imental indications from EELS |35| show no evidence for 
Al-Nb bonding, and suggest rather the existence of the 
O-terminated interface. According to our analysis this 
could only be marginally in thermodynamic equilibrium 
if the oxygen pressure is being 'buffered' by Nb/NbO and 
lies close to Pq" 1 * , which does not seem unreasonable. 

We have not included in our comparison interfaces 
with a different macroscopic orientation such as the 
Nb(110)/Al 2 O 3 (0001) interface @. Although this in- 
terface is believed to be thermodynamically more sta- 
ble than the Nb(lll)/Al 2 3 (0001) interface, the kinetic 
barrier to changing the macroscopic orientation is pre- 
sumably much greater than the barriers to changing the 
local interface structure. 

Although the formalism has been developed for de- 
scribing metal-oxide bonding, there are obvious applica- 
tions to systems in which water or other substances may 
contaminate surfaces or interfaces. The comparison of 
the energetics of interfaces with differing amounts of seg- 
regation follows the same lines. The application of the 
present formalism using a thermodynamic cycle to avoid 
the most difficult calculations may be fruitful in other 
situations in the field of interface chemistry. 
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TABLE I. Thermodynamic data used for calculating the 
dependence of surface and interfacial energies on oxygen par- 
tial pressure. The pressure represents the dissociation pres- 
sure of the oxide at 1500K. P° is one standard atmosphere. 





Q-AI2O3 


NbO 


AG" (kJ/mol) S 
log 10 (Po 2 /P°) 


-391.9 
-36.8 


-378.6 
-26.4 



TABLE II. In-plane relaxation of O and Al terminated in- 
terfaces and surfaces of alumina, a and Sr correspond to 
the angle of rotation and bond length increase of O triangles; 
surface and subsurface layers are indicated by the subscripts. 





N(b)/A/(0) 


N(b)/A(A1) 


A(O) 


A(A1) 


N(m)/A(0) 


Ql 


4.1 


4.2 


4.1 


3.1 


8.4 


OL2 


0.1 


0.2 


0.5 


0.2 


0.4 


Sri 


4.2 


4.3 


4.3 


3.2 


8.4 


5r2 


0.2 


0.2 


0.2 


0.5 


0.6 



TABLE III. W sep (in J/m 2 ) for both unrelaxed and relaxed 
structures. 



Interface 


Cleavage plane 


Unrel 


Rel 


N(b)/N(b) 


...Nb-Nb-Nb I Nb-Nb-Nb... 


4.9 


4.2 


A(A1)/A(A1) 


...Al-O-Al I Al-O-Al... 


7.0 


3.9 


A(A1)/A(0) 


...Al-Al-O I Al-Al-O... 


13.3 


12.7 


N(m)/A(0) 


Nb l O-Al-Al... 


10.9 


10.8 


N(b)/A(0) 


...Nb-Nb-Nb ; O- Al-Al-O... 


9.3 


9.8 


N(b)/A(A1) 


...Nb-Nb-Nb ; Al-O-Al-Al... 


4.2 


2.7 


N(b)/Nb-A(A1) 


...Nb-Nb-Nb I Nb-O-Al-Al ... 


4.0 


3.8 


N(b)/Al-A(A1) 


...Nb-Nb-Nb I Al-Al-O... 


4.4 


2.8 


N(b)-0/Al-A(Al) 


...Nb-Nb-0 I Al-Al-O... 


7.3 


4.9 



TABLE IV. Mulliken charges for atoms near the interface 
as a function of the excess of O. The notation for atomic 
planes is as in Fig. 1. 



r -s 


1.5 





-1.5 


Nb2 


0.27 


0.37 


-0.13 


Nbl/All 


0.77 


0.73 


0.36 


01 


-0.86 


-0.93 


missing 


A12 


1.49 


1.50 


0.15 


A13 


1.52 


1.51 


0.95 


02 


-1.00 


-1.00 


-0.99 
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FIG. 1. Side view of the Nb(lll)/Al 2 O 3 (0001) interface, showing labelling of the layers. 
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FIG. 2. Plane view of the (0001) surface of neutral alumina showing the lateral relaxation within the topmost O plane. The 
rotation and expansion of the O triangle below the surface Al atom is indicated by arrows. 
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FIG. 3. (a) Surface energies as a function of oxygen partial pressure at 1500 K. (b) Works of separation W sep (see Table III), 
(c) Interfacial energies, (d) Work of adhesion W a d, obtained by subtracting equilibrium interfacial energy from equilibrium 
surface energies. The region to the left between the vertical lines corresponds to the possible equilibrium states of the interface. 
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